Anomalous Josephson current via Majorana bound states in topological insulators 



P. A. loselevich^'^ and M. V. Feigel'man^-^ 
^L. D. Landau Institute for Theoretical Physics, Kosygin str.2, Moscow 119334, Russia and 
'^Moscow Institute of Physics and Technology, Moscow 141 700, Russia 

We propose a setup involving Majorana bound states (MBS) hosted by a vortex on a supercon- 
ducting surface of a 3D Topological Insulator (TI). We consider a narrow channel drilled across a 
TI slab with both sides covered by s-wave superconductor. In the presence of a vortex pinned to 
such a channel, it acts as a ballistic nanowire connecting the superconducting surfaces, with a pair 
of MBS localized in it. The energies of the MBS possess a 47r-periodic dependence on the super- 
conductive phase difference (p between the surfaces. It results in the appearence of an anomalous 
term in the current-phase relation, laif) for the supercurrent flowing along the channel between 
the superconductive surfaces. We have calculated the shape of the 47r-periodic function lai^p), as 
well as the dependence of its amplitude on temperature and system parameters. 



Since Majorana bound states (MBS) were predicted to 
exist in solid state systems, a number of different sug- 
gestions on how to detect them has been made. MBS 
are predicted to exist in systems characterized by both 
strong spin-orbit coupling and superconductivity. Exam- 
ples include the surface of a topological insulator (TI) 
covered by s-wave superconductor with vortices pj] or 
electrostatic defects [2|, as well as ordinary semiconduc- 
tor nanowires with spin-orbit coupling and proximity- 
induced superconductivity in the presence of a suffi- 
ciently strong Zeeman field Q. While there are numer- 
ous suggestions on detection of MBS [J, |3-[3 1 '^o exper- 
imental success has been reported yet. One more way 
to track down MBS in a superconducting proximity sys- 
tem is to observe an anomalous 47r-periodic phase depen- 
dence of a supercurrent. It was shown by A.Kitaev [8|, 
that a specific 1-D fermionic chain hosting a pair of MBS 
switches its ground state fermionic parity under an adi- 
abatic change of the superconducting phase difference Lp 
by 2-K. Thus, if the chain conserves fermionic parity, its 
behavior has to be 47r-periodic. 

In the present letter we propose and study a system 
with rather simple geometry, based on a TI sample cov- 
ered by a superconducting film, and calculate the anoma- 
lous 47r-periodic supercurrent it carries. More specifi- 
cally, we consider a flat thin slab of strong TI with both 
its surfaces covered by an s-wave superconducting film. 
The slab, together with the superconductive films, is 
pierced by a cylindric hole of radius i?, see Fig.l. A single 
superconductive vortex with a flux quantum $o = hc/2e 
is pinned to that hole. The two superconducting sur- 
faces are connected far away from the hole, forming an 
SNS circuit with the hole acting as a normal region. Be- 
low we demonstrate that MBS should be present in this 
setup, and calculate the anomalous 47r-periodic compo- 
nent la (</?) of the current flowing along the cylindric hole 
between the superconductive surfaces. Related issues 
were considered in a recent preprint Q in terms of the 
Hopf invariant and its relation to fermionic parity. Note 
that in the SNS-like setup we consider, a Hopf invariant 
cannot be defined since a part of the TI surface is not 




Figure 1: The system. A layer of TI has both surfaces covered 
by superconductor. A hole in the layer hosting a vortex forms 
an SNS-junction between the surfaces. The superconducting 
surfaces close away from the hole, completing an SNS-circuit 
with supercurrent flowing through the hole. 



gapped. 

In the simplest case of strong TI, realized in 8128 e^^ 
and Bi2Te3, surface electrons are described by a single 
Dirac cone with the Hamiltonian H = Vfa-p — Ef, where 
a is the spin operator and Ef is the Fermi energy of the 
surface states of the TI. Below we consider the semiclas- 
sical limit A ^ Ef and put h — 1. For the upper and 
lower surfaces covered by s-wave superconductor, a pair- 
ing term Aip'^ip'^ + h.c. arises due to proximity effect [l|: 

H = [vfCT ■ p - Ef )T^ + A(r)(rj; cos(p(r) -I- sin (^(r))(l) 

The Pauli matrices r act in the Nambu-Gor'kov space, 
K denotes complex conjugation. For the Hamiltonians 
Hu,i acting on upper and lower surfaces, one should re- 
place A(r) and (/j(r) in Eq.([T]) with A„^;(r) and </5u,;(r) 
correspondingly. The operator H anticommutes with the 
electron-hole conjugation operator S = ayTyK; we em- 
phasize that this property holds for any Bogolyubov-De 
Gennes Hamiltonian with a most general form of single- 
electron spectrum, including terms that break time- 
inversion symmetry. The symmetry {H, S} — di- 
vides the eigenstates of the Hamiltonian ([1]) into a set of 
conjugate pairs ijjE,i^-E = '^''Pe with opposite energies 
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{E, —E) and, possibly, a number of self-conjugate states 
with zero energy. In the basis of self-conjugated states all 
matrix elements of ([T]) are purely imaginary. The Her- 
mitian matrix H is antisymmetric in this basis, therefore 
the corresponding Pfafhan Pf H can be defined. It is 
easy to see that generically the existence of Pf H pro- 
tects twofold degenerate zero-energy levels of H against 
splitting under adiabatic variations of H. Indeed, con- 
sider a Hamiltonian with a parametric dependence H{ip) 
such that it possesses a pair of zero eigenstates at some 
ip = ipQ. Then the determinant of the matrix H has a 
double zero at ipo, i.e. (Pf H)"^ = det H{ip) (x {ip — ipoY , 
thus Pf H{ip) (X if — ipQ. Any small perturbation of 
the Hamiltonian (assuming it obeys the symmetry con- 
dition {H, S} = 0) can only shift the value of the cross- 
ing point (fioj but its very existence is robust. Appar- 
ently, the above arguments are in contradiction with 
well-known properties of Andreev levels in a usual super- 
conductive quantum point, with energies given jl4| by 

E± = ±AyJl — Tsin^ ^, so that an arbitrary small re- 
flection probability r — 1 — T leads to the splitting of zero 
levels present at (/s tt in the ballistic case T — 1- The 
origin of this contradiction lies in the spin degeneracy 
present in usual systems without spin-orbit interaction: 
in the presence of such a degeneracy the Pfaffian of the 
ballistic contact has a double zero, Pf H{(p) ^ [ip — tt)^, 
which is not robust to weak perturbations. Below we con- 
sider the generic case of strong spin-orbit coupling and 
thus no spin degeneracy. 

The sign of Pf H changes simulteneously with the 
fermionic parity Fq of the global ground-state of the sys- 
tem Indeed, the two eigenstates |e), |o) of the total 
Hamiltonian, which become degenerate when E^.o — 
al (p = (pQ, have fermionic numbers differing by 1. As 
the phase ip passes ipo, the roles of the ground state and 
the lowest excited state are interchanged. If the actual 
fermionic parity F of the system cannot change due to 
conservation laws (which we will assume to be the case), 
we come to the following conclusion: each time a pair 
of Majorana levels crosses E' = 0, the ground-state jgo) 
is transformed to the lowest excited state |eo), and vice 
versa, |eo) — ^ \9o)- We argue now that while the phase 
Lp changes on the (0, 27r) interval, an odd number of such 
crossings occurs, i.e. after a 2tt phase rotation our system 
does not return to its original state. 

Consider first the system shown in Fig. [1] without the 
cylindric channel but with two point vortices present in 
both superconductive films, on the upper and lower sur- 
faces. Each of them hosts a single MBS Due to 
finite thickness L of the slab, these two MBS xi,2 are hy- 
bridized into a single complex fermion V' = Xi +*X2- The 
energy eo{ip) of this fermionic mode is proportional, in 
general, to the amplitude of MBS tunnelling i ^ ^-^I^ti 
between the surfaces. However, for Lp = -k the tunnelling 
amplitude vanishes due to destructive interference (see 



Supplement 1), thus a single level crossing at -B = 
occurs as Lp varies on the (0, 27r) period. Let us now 
open the cylinder channel across the slab. It results in 
a drastic increase of hybridization between upper and 
lower superconductive surfaces, and in the appearence 
of an even (due to Kramers degeneracy) number of con- 
ductive modes. Since these additional modes appear in 
pairs only, the transformation — > eo occurs on each 
27r-period of Lp variation. The above arguments prove 
the existence of the anomalous component Ia{^) of the 
Josephson current which is odd under 27r-shift. Below we 
calculate its magnitude and temperature dependence. 

To find the current flowing along the hole channel we 
first calculate the sub-gap spectrum of the contact. We 
assume that the tube's radius R and length L (the latter 
coincides with the thickness of our TI slab) satisfy the 
conditions 

«i?,L«eo,6c (2) 

where ^0 = ^f/^^ Pf is the Fermi momentum of sur- 
face electrons of the TI, and ^sc = ■\/l?/2A <C ^0 is 
the actual " dirty-limit" coherence length in the supercon- 
ductive film with diffusion coefficient D. Inequlities ([2]) 
mean, in particular, that we consider a short SNS-contact 
with many transverse channels, Nch ~ PfR- To find the 
sub-gap spectrum we solve Bogolyubov-de Gennes equa- 
tions on both surfaces of the TI in the presence of the 
induced gap and vortices, and match obtained solutions 
with eigenmodes living on the inner cylindric surface of 
the channel. 

We use cylindric coordinates r, 0, z with the z-axis co- 
inciding with the tubes axis. The Hamiltonian ([1]) in the 
presence of an Abrikosov vortex is: 

H = vf<j- (pr, - V0/2) - Eft, + A(r)T, (3) 

The magnetic screening length is very long for thin films, 
and the flux of the vortex is distributed over a radius 
much greater than both R and ^0, thus we may neglect 
the vector-potential term in ([3]) . For a fixed angular mo- 
mentum 1/ we get = Qi^o-icr^0/2 I "('"M with radial 

\v{r)) 

wave functions u[r) and v{r) given by 

with H„i standing for Vf [a^ [pr - ^) + ctvt) ~Ef. The 
profile of the gap function near the vortex center is 
A(r) ?a rA/y/r^ + 2^^^. Equation (|H) can be solved for 
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e ^ A In the first order in e/A it yields: 



(5) 
(6) 



0(r) 



,,(1.2) 



K(r) = / Aip)dp (7) 
Jo 

^'^^-.Ki,, (8) 



(r) 



(9) 



H^^''^^ in ([9]) denote Hankel functions. Wlien iy,e — 0, 

Due to tlie fac- 
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solutions dS]),® become exact 
tor in Eqs.lIS]),®, the wavefunction is localized in the 
vicinity of the tube within a length ^o- Next we specify 
electronic eigenmodes in the channel. The Hamiltonian 
for a cylindrical surface has the form 12 1 



(JzPz 



2R 



(10) 



with ae = —<Jx sin + cos 9 and ar — Ux cos 9+ ay sin ( 
The eigenfunctions of (|T0| are 



COS ■ 



(11) 



where a = arctan-^. The energy spectrum for positive 

energies (counted from Ep) is e^ p = \Jp'^ + The 
angular momentum is half-integer, jj, = ^ + n, due to the 
Berry phase originating from the rotation of the spin. 
We neglect the effect of the small magnetic flux ^ 
<i>o inside the cylinder. Note that the possibility to find 
conductive channels with no magnetic flux inside is due 
to the semiclassical condition pfR^ 1; on the contrary, 
in the ultraquantum limit PfR < 1 the presence of a $o 
flux would be necessary to make the hole conductive, see 
Supplement to Ref. [13| . 

Now we have to match wavefunctions ([S]),® on the 
upper and lower surfaces with wavefunctions (jlip (and 
their analog for the hole component with e < 0) on the 
boundaries between the cylinder and fiat surfaces. The 
proper matching condition for a sharp edge reads: 



'Si — exp 



-^-(cr-ne) 



(12) 



Here is the unit vector in the direction of the edge and 
9 is the rotation angle (equal to ^ in our case). ^1^2 are 
4-component vector wavefunctions on both sides of the 
edge. The operator in Eq. (fT2l) rotates the spin by the 
angle 9 around the direction ne. Wc provide an explicit 
form of the matching equations and the derivation of the 
resulting energy spectrum in the Supplements 2,3. 
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Figure 2: The anomalous current lai^'f) is computed for 
PfR = 2 (dashed line) and pfR — 3 (blue line). Other pa- 
rameters are fixed as p/^ = 10, p/^sc = 5, PfL — 6, T = 
0.05A. The extrema of laif) occur whenever a pair of con- 
jugated levels crosses at e = 0. 



The resulting low-lying levels with e ^ A as function 
of the phase difference ip are given by (up to the neglected 
terms of the order of Se ^ vA/{pfRY): 



PfR 2p/i?2 



27rk) - 

-1/2- 



1 - 



[pfRY 



(13) 



where Ei « A and eo ~ A/2 for the case ^sc ^ £.0 we 
consider, < pfR and k is an integer. No propagating 
modes in the tube exist aX v > pfR, instead there is an 
exponentially small overlap between the bound states re- 
siding inside vortex cores on opposite surfaces. The ener- 
gies of the global eigenstates with \v\ > pfR are given by 

= t^^;^ log ^ up to exponentially weak (/s-dependent 
corrections. The result ((T5)) is not applicable in the re- 
gion \iy — pfR\ ^ {pfRy^^ where the crossover between 
hybridized and non-hybridized levels occurs. 

The supercurrent through a short SNS-contact in ther- 
mal equilibrium can be expressed [l3 | in terms of An- 
dreev levels: I{ip) = — 2e ^^-.^^.^g tanh (^^7) Here 
all eigenstates are taken into account, regardless of their 
parity. The total current scales with the number of con- 
ductive channels in the hole: I ~ eApfR/Ti. To re- 
veal any parity-related effects, we have to consider ther- 
modynamic ensembles with odd and even numbers of 
quasiparticles separately. A standard route is to in- 
troduce thermodynamic potentials ^odd/even, describ- 
ing odd and even numbers of quasiparticles correspond- 
ingly [l3|. Dividing the total current I[ip) into a sum 
Iniv) + laif) where /n(<p) is parity- independent, we ob- 
tain/a = (-l)^"e^(riodrf-f^ei;e«), seeSuppl.4. In terms 
of the spectrum of Andreev levels, it reads 

2e/ 1 dej 



U^) = (-1)^" 



(14) 



where / = ^ - tanh |y = jhy\, ■ fnon and both the prod- 
ucts and the sum are done over levels with e, > 0. Fac- 
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crossover temperature T2; however, we do not expect 
their effect to be drastic. Finally, the usual parity- effect 
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Figure 3: Temperature dependence of /a(</' = for the same 
sets of parameters as in Fig[51 Blue and dasfied lines corre- 
spond to pfR = 3 and pfR — 2. 



tors fhyb and /„o„ correspond to the hybridized (current- 
carrying) and non-hybridized levels. The physical mean- 
ing of la is the current difference between an odd and an 
even state of the system. Eqs. (|13|14|) constitute our ma- 
jor quantitative result. Two examples of the /a (95) depen- 
dence computed using Eqs. (jf 3lf 41) are presented in Fig. [21 
Ia{'p) experiences » 2p/i? oscillations in the (— tt, -t-7r) 
interval and has opposite signs at ip = ±7r where the 
amplitude of lai^) is maximal. Note that approxima- 
tions used to derive Eg. dT^ may lead to a deformation 
of the /a(<p) dependence, leading to an inhomogeneous 
shift of its oscillating pattern by the amount ~ {pfR)~^. 
It should not affect, however, the maximal value I^'^^ of 
the anomalous current. The temperature dependence of 
jmax jg presented in Fig. [3] 

To analyze the temperature dependence of laif) we 
consider Eq. p^ in several limits. In the range of T much 
higher than the typical level spacing ei/pfR, the sub- 
product fhyb ~ TpfR/iei ^f-Qj. ^j^g derivation of this 
and the following formulae see Supplement 4). If, in ad- 
dition, T ;» £'0 = In then the subproduct /„o„ is 

also small: fnon — e^^^"^/^*, where 5 — {A/pf(o) In is 
the typical level spacing within a single vortex core. Thus 
la decays exponentially with temperature at T > Ti — 
min(^, ^ In At lower temperatures T < Ti we 

find fhyb = tanh(ei/2T), where £i{(p) is the lowest hy- 
bridized level, therefore 



eeo • fno 



cosh^ - sinh^ ^ 



fnon — GXp 



2T non '^^^^^^ 2T 

4e-£o/T 



1 - e-VT 



(15) 
(16) 



Simple analysis of Eas. (|15ll6p leads to the second char- 
acteristic temperature T2 = Eq/ ln{pfS^o). Depending on 
parameters, T2 may be both higher and lower than Ti. 
In addition, we mention the existence of non-hybridized 
subgap states localized in superconductive films near vor- 
tex cores, which may lead to some suppression of the 



temperature T3 — A/\ii{iyVA), see Ref. [15|, puts an 
additional restriction for the temperature region where 
an anomalous current could be observed. 

Summarizing the above analysis, we find that for the 
anomalous current la to be detectable, the following con- 
dition must be met: 



T < min 



A 2Ai?ln(eo/6c) 



PfR' ^oHPf^o) ' H'^VA) 



(17) 



The temperature dependence of the anomalous current 
is shown in Figl3] for a specific choice of parameters such 
that Ta < Ti. 

To conclude, we proposed a setup using strong topolog- 
ical insulator covered by superconductive films, which al- 
lows the detection of Majorana bound states through the 
measurement of an anomalous 47r-periodic component la 
of the Josephson current. The temperature dependence 
of the la amplitude is calculated, and the conditions for 
the proposed effect to be observed are found. 
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the RFBR grant # 10-02-00554 and by the RAS program 
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ONLINE SUPPLMENTARY MATERIAL 
1. Degeneracy of MBS at ip — n 

Beside of the S-symmetry inherent to Bogolyubov-de 
Gennes equation, our model has a symmetry that is the 
composition [/ of rotation i? and time-reversal 0. Indeed, 
by rotating our system around an axis lying in the z — 0- 
plane (the plane lying between the surfaces of our slab), 
we arrive at the same system as by conjugating A, pos- 
sibly with a phase shift. Selecting a gauge, in which the 
superconducting phase is given by argA(0, z ~ L/2) = 
9 — arg A(6', z — L/2) — ip, we find that U commutes with 
the BdG Hamiltonian if the axis of rotation forms a (p/2 
angle with the x-axis. The part of R that acts on spin 
variables equals — i[ax cos ip/2 + ay simp /2], while O — 
icTyK, where K stands for complex conjugation. Hence, 
U ~ sin ip/2 + iGz cos <y9/2. The first term commutes with 
S = GyTyK, while the second term does not, in general. 
We see that U commutes with S only if = tt. The 
commutation relations [H, U] = 0; [C/, S] =0; {if, S} = 
lead to the degeneracy of the E — Q levels of our system 

at ly9 = TT. 



2. The matching equations 

Consider two surfaces - 1 and 2 - joining on a line 
at the angle of ^o- Let ^'1,2 be the wave function on 
the edge of surface 1 and 2 respectively. We replace the 
sharp boundary by a smooth cylindric transition with a 
radius r. We use the natural cylindric coordinates of this 
cylindric sector with = at its boundary with surface 
1. Next we decompose ^'1^2 into the set of functions (11): 



^2 =^ap,^e'''''-*'^^''°/2^'p,, 



(18) 
(19) 



with 



In the limit r ^ Q we have to 



choose /X = since \J + ^ = e^,^ = const. This allows 

us to write '^2 = e^"^^^°/^5'i. The basis-independent 
form of this equation is given by Eq.(12). 

3. The spectrum equations 

The set of equations defining our systems spectrum 
consists of 8 equations - two sets of 4 matching equations 
corresponding to the two ends of the cylindric hole. To 
derive these equations, let we first write down the wave 



functions for some fixed momentum v; we will do it in 
linear approximation in small parameter <C 1. The 
fixed angular momentum v on the surfaces corresponds to 
v± = J/ ± 1/2 in the channel: electron waves have /i = i/_ , 
while hole waves have fj, = After we have chosen the 
proper angular momenta, we only need to match wave 
functions at some fixed angle, say 9 — 0. The w.f. in the 
cylinder is a superposition of four waves, representing 
electrons and holes propagating up and down the tube: 



— g«(P/+9)zcosQ_ 



-i{pf+q)zc 



a± 



-i(pj—q)z cos 

v±l/2 



cos Q!_/2 

i sinQ!-_/2 
-zsina_/2 

C0SQ!_/2 

cosq;+/2 
isv[\a^/2 

-i sina^ /2 
cos a+ /2 



arcsm ■ 



PfR 



(20) 
(21) 
(22) 
(23) 
(24) 



Next we write out the surface wave functions. At r ^ ^ 
the solution ^ transforms into 



^u-l,2i^+)-^^^''^"^u-l,2i^+) (25) 
-#o,„(2) 



V = -ic,e-w^l,i^{r,) + C2e-*^°<^,/2(r_) (26) 



with r± 



r{l ± q/pf) and 



h = 0(0) ~ g?o/2 + 



^^j^ log ^ in the dirty limit Co- Using the asymp- 

totics of Hankel functions with large arguments, we find 
the w.f. on the upper surface at r = R, z = L,9 = 0, 
within ^ 1/pfR accuracy (we neglect terms ~ {pfR)^"^): 



^upiR) = Ci 



+C2 



i(j)o + ip] 



/ exp[i/3 — 2i(j)o 
i exp[— 2i0o 
1 

y iexp[—i/3] 

( exp[-i/3 + 2i(f>Q + ip\\ 
i exp[2i(/)o 
1 

i exp[i/ 



(27) 



with some coefficients Ci,C2 and cos/3 = v/pfR. This 
formula is correct unless \v~pfR\ < [pjR)'^/'^. '^down{R) 
is obtained by introducing two new coefficients Qi,Q2 
instead of Ci, C2 and putting (p = 0. 

Now we can construct the matching equations accord- 
ing to the matching rule (12). Acting on '^cyi{z = L) and 
"^cyiiz = 0) by ^(1 - iay) and ^(1 + icTy) respectively 
we finally get 
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e 

1 



A, 






h 



= Cl 



exp[il3 — 2i(j)Q + ip] 
P 
1 



z exp[— 

/ exp[i/3 - 2i(;/)o] 
^ I z exp[— 2j<j!>o] 



i exp[— j/3] 



i exp[j/3] 



exp[-i/3 + 2i0o] 
j exp[2i(^o] 



i exp[i/3] 



Since a± — 7r/2— /3 within our approximation, the spinors 
with coefficients A-^ and Ci in the first equation of the 
system (|28p are collinear. The same is true for the spinors 
with coefficients and C2 in the same equation. A 
similar statement is true for each of the other Egs. ipSl) . 
Consequently, the system (1281) splits into two simple sub- 
systems: 



A-^ = 

Bi = 

^f- = 

Bi = 



Cl exp[i/3 — 2icj>o + i(p — i{pf + q)L cos a-] (29) 
Cl exp[— i(p/ — q)L cos a+] 
Q2 exp[-j7r/2 + 2i(l>o] 
Q2 exp[-j7r/2 + j/3] 



and 



where the upper/lower signs refer to the systems 
and pop correspondingly. Since we assume L <C $0 and 
i? ^ we neglect the ~ qL term in the r.h.s of (|3T]) and 
the j'-dependent term in 0o i arriving at the final formula 
(13). 



4. The anomalous current 

The thermodynamic potential of a system with a fixed 
parity can be written as 



Ai = 
Ai = 



' exp[j7r/2 + 2i</>o + i(p + i{pf + g)LcosQ_] (30) 
C2 exp[j7r/2 + if} + i{pf — q)L cos q+] 
Qi exp[— 2j(/!>o + ijS] 



These equations give the following set of the two spectral 
equations 



400 = -2a ± ((/? - TT - 2'Kk) + 
-{pf — q)L cos a+ + (l)L cos a_ 



(31) 



^odd/even — T hi 



-(32) 



with (3 — T ^. The difference 5Vl = ftodd ~ ^even equals 



Sn = -T In 



1-/ 
1 + / 



Hence 



with / = J|tanh^ (33) 



" ^ " 1 - " 1 - /2 4^ 2cosh2 4|-tanh^ 



1 _ « 



j2 sinh/3ei 



(34) 



r 



It is useful to divide / into factors corresponding to 
the hybridized and non-hybridized parts of the discrete 
spectrum. 



We neglect the exponentially small dependence of non- 
hybridized energies on (/?, so that 



^^fhybfno 

1 _ f2 72 

•I hybJ non 



-J 



J 



, sinh(e./T) 



(36) 



/ = Y[ <^anh 1^ Jl tanh ^ 



W\>PfR 



2T 



fhybfnon (35) 



First let us analyze fhyb- The typical level spacing 
is Wo If r > loq, we can write = 



exp lntanh(^^)(i:^] = exp 



4wo 



For small L/R 



7 



and ip = Tr this gives exp — ^ 2A^^ opposite case 

T ^ cjo wc can write fhyb = tanh^, where is the 
lowest hybridized energy level. 

Next we consider factor J. At T ^ ojq we use 

sinh(e/T) ~ e/T at small energies and get J < 

EI-? ^ = 2:^log^. At T < Wo the amplitude of 
J is dominated by the lowest hybridized level ei and we 



get |J| = 2sintjT - 

Finally, we analyze fnon- The spectrum of the non- 
hybridized de Gennes states is described by Enon = 
^ln(^) \pfR + n] = Eo + n6 with n = 0,1,2.... 

There are two copies of these series one for the up- 
per and one for the lower surface, li {T — Eq) ^ 6 we 
can rewrite 



fnon — exp 



PfR 



1 1 '^'^ , 

in tann — an 
2T 



exp 



26 



In tanh xdx 



(37) 



The second term can be neglected, ifT ^ Eq. At low temperatures T <^ Eq we have 



exp 



exp 



. _£0 \ ^ _riS_ 

-4e T > e T 



n=0 



exp 



e T 
1 — e T 



(38) 



The low-temperature expressions derived for fhyb, fnon, J lead to formulae (15,16). 



